arXiv:1504.07476vl [gr-qc] 28 Apr 2015 


Warm intermediate inflation in the Randall-Sundrum II model in 


the light of Planck 2015 and BICEP2 results: A general 

dissipative coefficient 

Ramon HerreraQ 

Instituto de Fisica, Pontificia Universidad Catdlica de Valparaiso, 
Avenida Brasil 2950, Casilla 4059, Valparaiso, Chile. 

Nelson VidelaQ 

Departamento de Fisica, FCFM, Universidad de Chile, 

Blanco Encalada 2008, Santiago, Chile 

Marco OlivareJil 

Facultad de Ingenieria, Universidad Diego Portales, 

Avenida Ejercito Libertador 44^ Casilla 298-V, Santiago, Chile. 

(Dated: April 29, 2015) 

Abstract 


A warm inflationary Universe in the Randall-Sundrum II model during intermediate inflation 
is studied. For this purpose, we consider a general form for the dissipative coefficient T{T,4>) = 
C(j) ; and also analyze this inflationary model in the weak and strong dissipative regimes. We 
study the evolution of the Universe under the slow-roll approximation and find solutions to the full 
effective Friedmann equation in the brane-world framework. In order to constrain the parameters 
in our model, we consider the recent data from the BICEP2-Planck 2015 data together with the 
necessary condition for warm inflation T > H, and also the condition from the weak (or strong) 
dissipative regime. 
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I. INTRODUCTION 


In modern cosmology, our notions concerning thejplwsics of the very early Universe have 


introduced a new element, the inflationary period 


which provides an attractive ap¬ 


proach for solving some problems of the standard Big-Bang model, like the flatness, horizon, 
monopoles, among other. However, the essential feature of inflation is that can generate a 
novel mechanism to explain the Large-Scale Structure (LSS) of the Universe j7Hll[| and 
provide a causal interpretation of the origin of the anisotropies observed in the Cosmic 


Microwave Background (CMB) radiation 12l-ll8| 


By other hand, with respect to the dynamical mechanisms of inflation, the warm inflation 
scenario, as opposed to the standard cold inflation, has the attractive feature that it avoids 
the reheating period at the end of the accelerated expansion, because of the decay of the 
inflaton held into radiation and particles during the slow-roll phase 19|. During the warm 
inhation scenario the dissipative effects are important, since radiation production occurs 
concurrently together with the inhationary expansion. The dissipating ehects arise from a 
friction term which describes the process of the scalar held decaying into a thermal bath. 
Also, during the scenario of warm inhation the thermal huctuations arises from the inha¬ 
tionary epoch may play a fundamental function in producing the initial huctuations indis¬ 
pensable for the Large-Scale Structure (LSS) formation [20, l2l|. An indispensable condition 
for warm inhation scenario to occur is the presence of a radiation held with temperature 
T > H during the inhationary expansion of the Universe. Since the thermal and quantum 


huctuations are proportional to T and H, respectively 19l-l21|. so when T > H, the thermal 


huctuations of the inhaton held predominates over the quantum huctuations. Also, inhation 
ends when the Universe heats up to become radiation dominated, and then the Universe 
smoothly enters to the radiation dominated era, without the need of a reheating scenario 
For a review of warm inhation, see e.g. Ref. 221. 

On the other hand, cosmological implications of string/M-theory have currently attracted 
a great deal of attention; specihcally some were concerned with brane-antibrane conhgura- 
tions such as space-like branes j^. In this conhguration, the standard model of particles 
is conhned to the brane, while gravitation propagates into the bulk space-time. Here, the 


ehect of extra dimensions induces extra terms in the Friedmann equation 


24 


26| . In partic¬ 


ular, the cosmological Randall-Sundrum H model (RS H) has received great attention in the 
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last years [27|. Randall and Snndrnm snggested two similar but phenomenologically different 
brane world scenarios. In this form, there are two versions of the Randall-Snndrnm model, 

ge„e.all. .en.oned . Ra.dall-Sund.u. . (RSl)S aad Raadall-Sund.a. 1. (RS2,0. 

Somewhat confusingly, the RSI model includes two branes, while the RS2 model only con¬ 
tains a single brane. 

In the RSI model, we have two 3-branes separated by a region of five-dimensional anti-de 
Sitter spacetime. Here, the fifth coordinate is compactified on S'^/Z 2 , and the branes have 
equal; positive and negative tension branes are on the two specified points. In this model 
of brane the matter fields are confined on the two brane and the gravity propagates in the 
5-dimensional bulk. However, one important point of discussion in this brane model, is the 
mechanism that works to select the neces sary separation distance between the two branes 


called radius stabilization 


291 (see also Ref.|30|). Similarly, the stable 


ization mechanism plays 


3ll. 


a crucial role in the recuperation of 4-dimensional Einstein gravity 

In contrast, the RS2 modelj^ contains a single, positive tension brane and a non-compact 
extra dimension is infinite in extent and the radius stabilization is not present in this brane 
model. By other hand, in the RS2 model, the observable universe is a four-dimensional 
time-like hypersurface. The five-dimensional energy-momentum tensor can be spiffed in a 
regular and a distributional part. The regular part describes the non-standard matter fields 
in five dimension and the distributional term containing the brane tension and the standard 
matter fields on the brane. Except gravitation, all standard model interactions and matter 
fields are confined to the brane. 

These alternatives to Einstein’s General Relativity cosmology belongs to the so called 
brane-world cosm olog ical models. For a comprehensible review of brane-world cosmology, 


see e.g. Refs. 


32 


35 


In the observational aspect, nowadays there is strong evidence 


that the very early Universe could have experimented an inflationary expansion period, as 
was pointed out in the begging of the introduction. An important feature of the inflationary 
period, is that is located in a period of cosmological evolution where the effects predicted 
by string/M-theory are relevant. For this reason, over the last decade, there has been great 
interest in the construction of inflationary models inspired in these theories. In the following 
we will concentrate only in the RS2 model, which forms the basis for the rest of this work. 

As regards exact solutions, one the most interesting in the inflationary Universe can be 
found by using an exponential potential, which is often called a power-law inflation, since 
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the scale factor has a power-law-type evolution, here the scale factor is given by a{t) ~ t^, 


in which p > 1 


37l |. Also, an exact solution can be obtained in the de Sitter inflationary 


Ij. Moreover, exact 


Universe, where a constant effective potential is considered; see Ref. 
solutions can also be found for the scenario of intermediate inflation, for which the scale 
factor a(t) evolves as 

a{t) = exp[At-^], ( 1 ) 

where A and / are two constants; A > 0 and 0 < / < 1 


It is well know that the 


expansion rate in this model is slower than de Sitter inflation, but faster than power-law 
inflation; this is why it is known as ’’intermediate”. This inflationary model was originally 
developed as an exact solution, however, intermediate inflation may be described from the 
slow-roll analysis. Under the slow-roll approximation the scalar spectral index n* ~ 1, and 
for the specihc value of / given by / = 2/3 the spectral index is = 1 (Harrizon-Zel’dovich 
spectrum) 39|]. Also, an important observational quantity in this model, is that the tensor- 
to-scalar ratio r, which is signihcantly r 7 ^ 0 40|, 1^. Other motivation to consider this type 
of expansion for the scale factor comes from string/M-theory, which appears to be relevant 
in the low-energy string effective action l43(] (see also. Refs. 44l-l47l| ). These theories can 
be utilized to resolve the initial singularity and describe the present acceleration in the 
universe, among others 48(]. Also, the approach of consider the warm intermediate inflation 
has already been studied in the literature, ^ in the context of other frameworks. 

In this way, the goal of this paper is to analyze the possibility that a higher dimen¬ 
sional scenario, in particular the RS 11 brane-world model, can describe the dynamics of 
the Universe in its very early epochs. We propose this possibility in the context of warm 
inflation scenario for a Universe evolving according to the intermediate scale factor, and how 
a generalized form of dissipative coefficient r(T, 0) oc inffuences the dynamics of 

our model. In order to study our brane-warm intermediate model we will consider the full 
effective Friedmann equation, and not the lower energy limit or the high energy limit for 
the effective Friedmann equation. Also, we will study the cosmological perturbations, which 
are written in terms of several parameters. Here, the parameters are constrained from the 


B1CEP2 data 


50| together with the Planck satellite 17| and Planck 2015 18|. By other hand. 


it is well known that the B1CEP2 experiment data has important theoretical signihcance on 
the amplitude of primordial gravitational waves produced during inflation. In this form, the 
tensor-to-scalar ratio r from the B1CEP2 data, has been found at more than 5a conhdence 
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level (C.L.) in which the ratio r = 0 . 2 ]'^g;Q 5 at 68 % C.L., and with foreground subtracted 
r = 0.16^'^n n® 50|]. Nevertheless, this value of the tensor-to-scalar ratio has become less 


^- 0.05 

transparent , see e.g. Ref. 


51| . Of this way, a detailed analysis from Planck and BICEP2 


data would be necessary for a definitive answer to the diffuse Galactic dust contamination. 
Recently, the Planck collaboration published new data of enhanced precision of the CMB 
anisotropies}]^ . Here, the Planck full mission data improved the upper bound on the tensor 
to scalar ratio ro .002 < 0.11(95%C'.L.), and this upper bound for the ratio r is similar to 
obtained from Refs. [l3, 0, in which r < 0.12(95%G.L.). 

The outline of the paper is as follows: The next section presents a short review of the 
effective Friedmann equation for the Randall-Sundrum type II scenario. In the section III 
we study the dynamics of warm inflation in this brane-world model, in the weak and strong 
dissipative regimes; specihcally, we obtain explicit expressions for the inflaton scalar held, 
dissipative coefficient, and effective scalar potential. Immediately, we compute the cosmo¬ 
logical perturbations in both dissipative regimes, obtaining expressions for observational 
quantities such as the scalar power spectrum, scalar spectral index, and the tensor-to-scalar 
ratio. At the end, section IV summarizes our results and exhibits our conclusions. We chose 
units so that c = h = 1. 


II. THE BRANE-INTERMEDIATE WARM INFLATION SCENARIO 


Followings, Refs. 25|, l26| we consider a hve-dimensional scenario, for which the modihed 
Friedmann equation for a spatially hat Friedmann-Robertson-Walker (FRW) metric, has the 
form 


= up 


1 + 


2 r. 


A4 

T 


( 2 ) 


where H = a/a represents the Hubble parameter, a denotes the scale factor, the constant 
K = SttG/S = Sn/Snip, and the energy density p corresponds to the energy density of 
the matter content conhned to the brane, A 4 is the ehective four-dimensional cosmological 
constant on the brane, and denotes the inhuence of the bulk gravitons on the brane, 
where ^ is an integration constant. The term r represents the brane tension and, considering 
the nucleosynthesis epoch, the value of the brane tension is constrained to be r > (IMeV)^ 

. However, a stronger limit for the value of the brane tension results from usual tests for 
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deviation from Newton‘s law, for which r >(10 TeV)"^, see Ref. 341. 

In the following, we will consider that the cosmological constant A 4 = 0, and once inflation 
begins, the term rapidly becomes unimportant. In this form, the effective Friedmann 
equation given by Eq.([2]) becomes 


= K, p 


1 + ^ 

2t. 


(3) 


On the other hand, during warm inflation the Universe is filled with a self-interacting 
scalar field with energy density p^j) together with a radiation field with energy density p-y. In 
this form, the total energy density p of the Universe can be written as p = p^p + p^. Here, the 
energy density pp and the pressure of the scalar field are given by p^p = 0^/2 -|- V (0) and 
P<l> = 0^/2 — U(0), respectively. The term U(0) represents the effective scalar potential. In 
the following, we will consider that the dots mean derivatives with respect to cosmic time. 

We will assume that the total energy density p is confined in the brane, and then the 
continuity equation for the total energy density becomes p + 3 H {p + P) = 0. In this 
way, following Ref. , the dynamical equations for and p^ during warm inflation can be 
written as 

Pp + 3 H (pp Pp) = —r 0^, (4) 

and 


p^ + AH p^ = r0^. 


(5) 


Here, F > 0 represents the dissipative coefficient and describes the process of scalar field 

le inflationary expansion 1^. In the context of brane- 


52l | was studied a high energy scenario during the strong 


decaying into radiation during t 
warm inflationary model, in Ref. 
dissipative regime. 

From quantum field theory, the dissipative coefficient F was computed in a supersymmet¬ 
ric model for a low-temperature scenario In particular, for a scalar field with multiplets 
of heavy and light fiel^ it is possible to obtain several expressions for the dissipative coef¬ 


ficient F, see e.g.. 
Following Refs, 
given by 




53-58 


57( |. we consider a general form for the dissipative coefficient F(T, 0) 




( 6 ) 
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where the constant is related with the dissipative microscopic dynamics and the expo¬ 
nent m is an integer. This expression for the dissipative coefficient inclndes different cases, 


depending of the valnes of m, see Refs. [^,1^. Concretely , for the special valne of m = 3, 
for which T = the parameter agrees with = 0.02h^A/y, where a generic 

snpersymmetric model with chirial snperfields <h, X and Yi, i = I, ...Afy is stndied [53, l59 |. 
For the special case m = 1, the dissipative coefficient is related with the high temperatnre 
snpersymmetry (SUSY) case. Finally, for the cases m = 0 and m = —1, the term F(T, 0) 
represents an exponentially decayin g pr opagator in the high temperatnre SUSY model and 


the non-SUSY case, respectively 




Considering that in the scenario of warm inflation the energy density predominates 
over the density then the Eq.([n]) becomes 






1 + 


02 + 2U(0) 


4 r 


(7) 


In the following, we will not stndy the effective Friedmann eqnation, given by Eq.(I7j), in 
the lower energy limit i.e., <C r or in the high energy limit i.e., S> r as onr starting 
point, instead we will consider the fnll effective Friedmann Eq.d?]). Here, we note that 
there are two ways of deriving the Friedmann’s eqnation from hve-dimensional Einstein’s 
eqnations. The hrst method is rather simple and considers only the bnlk eqnations. The 
second approach ntilizes the geometrical relationship between fonr-dimensional and five- 
dimensional qnantities. However, the Einstein eqnations and in particnlar the Friedmann’s 
eqnation in the bnlk inclnde different fnnctions (from FRW metric in 5-dimensional). In 
particnlar, these fnnctions are subjects to conditions (junction conditions) on the brane 
localized at p = 0 and symmetry (Z 2 -symmetry) when integrate over y. In this form, we 
could not obtain analytical solutions considering the full 5-dimensional equations of motion 
from 4-dimensional analytical solutions. In the following, we will obtain analytical solutions 
in four-dimensional of the full effective Friedmann’s equation only. 

In this form, solving the cuadratic equation ([7]) for p^ (where we take the solution for 
which < 0), and combining with Eqs. (jl]), results 


-H) 


1 + 


2H^ 


KT 


- 1/2 


^ 3k (1 + R) 

where the quantity R denotes the ratio between F and the Hubble parameter, dehned as 

F 


( 8 ) 


R = 


3H' 


(9) 
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In the following, we will consider that for the case of the weak or strong dissipation regime, 
the ratio R satishes R <1 or R > 1, respectively. 

Also, following Refs. 


19 


20| . we consider that during warm inflation the radiation pro¬ 
duction is quasi-stable, i.e., and p^ r0^. In this form, combining Eqs.(|5]) and 


), the energy density of the radiation held results 

r02 V{-H) 

P-y 


1 + 


2H^ 


K T 


- 1/2 


( 10 ) 


AH BkHII + R) 

On the other hand, we assume that the energy density of the radiation held is given by 
p^ = C^T^, where the constant = 7r^5f*/30 and the term represents the number of 
relativistic degrees of freedom. Combining the above expression for the energy density p.y 
and Eq. flTOll . we hnd that the temperature of the thermal bath T yields 


T = 


-H) 


-| 1/4 



2 m 


1 1 


KT 


6 K C^H (1 -|- R) 

From Eqs.dnj and ([H]), the ehective scalar potential results 


- 1/8 


( 11 ) 


V = T 


-1+t 1 + 


2 m 


+ 


H 


KT j 3«:(1 -|- R) 


1 + 


2 m 


2 \ - 1/2 


KT 


( 12 ) 


Here, we note that the ehective potential could be determined explicitly in terms of the 
scalar held 0, in the weak (or strong) dissipative regime. 

Now, combining Eqs. lfTT]) and (jb]) we get 

m/4 


_ 4 —m _ 

r 4 = c^<p 


1—m 


-H 


6k C-yH{l + R) 


1 + 


2 m 


2 \ -m/8 


KT 


(13) 


We note that the above expression determines the dissipation coefficient in the weak (or 
strong) dissipative regime in terms of the scalar held 0. 

In the following, we will study our brane-model for the general form of the dissipative 
coefficient, given by Eq.(|6]), during intermediate inhation, where the scale factor evolves 
according to Eq.(IT]). We will restrict ourselves to the cases R < 1 (weak regime) and R> 1 
(strong regime). Also, in the following, we will study the cases m = 3, m = l, m = 0 and 
m = — 1 corresponding to the dissipative coefficient E = r(T, 0). 
























A. The weak dissipative regime. 


Firstly, we consider that our brane-model evolves agreeing to the weak dissipative regime, 
i.e., R < 1 (or analogously F < 3H). In this from, combining Eqs. llT]) and ([H]), the solution 
for the standard scalar held 0 results 

B[t] 


- 00 = 


K 


(14) 


here the quantity (j){t = 0 ) = 0 o is an integration constant, than can be assumed as 
0(t = 0) = 00 = 0 (without loss of generality), and the constant K is specihed by 
K = . The function B[t], represents the incomplete Beta function 

Krf2(i-/) 1 : 


Af ITPF 


6 l| . given by 


B[t] = B 


2(A/)2 ’4(l-/)’4_ 

From Eqs. ffT]) and m, we hnd that the Hubble parameter H as function of the scalar 
held, results -ff(0) = A f {B~^[K 0]) ^ . Here, the function B~^[K 0] corresponds to the 

inverse of the incomplete Beta function B[t] 61|. 

The ehective potential as function of the scalar held 0, from Eq. fll2p and considering that 
V > 0^/2, results 


1 /( 0 ) ~ r 


■1 + 11 + 


2A‘^P{B-^[K(j)]) 


- 2 ( 1 -/) 


K T 


(15) 


Considering that i? < 1, from Eq. flT^ . the dissipative coefficient F in terms of the scalar 
held becomes 


r( 0 ) = c; 


4 

4 —m 


1 -/ 


_QkC^B-^[K (t)\_ 


4-m 4(l-m) 

(5 4-m 


1 + 


2/l2f (H-i[iC0]) 


- 2 ( 1 -/) \ 


K T 


, (16) 


for the case m 7 ^ 4. 


From the dehnition of the dimensionless slow-roll parameter e = —jp, then we hnd that 

£ = form, and considering that the condition for inhation to occur 

is determined by e <1 (or equivalently a > 0 ), then the scalar held during warm inhation 

1 r/i_f \ 

satishes the condition (j) > j^B ( -jj- j 

On the other hand, from Eqs.([T]) and flTT)) . the number of e-folds N between two diherent 
values of the scalar held, denoted 0 i and 02 , is 
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N = £ Hdt = A = A [{B-^[K (j)2]y - {B-^[K (j)^]y] . (17) 

Following Ref.|4]J, the inflationary phase begins at the earliest possible scenario. In this 
way, the scalar held (pi takes the value 


b = —B 
^ K 


f 


Af 


1 // 


(18) 


In the following, we will study the scalar and tensor perturbations for our brane-warm 
model during the weak dissipative regime. For the case of the scalar perturbation, it could 
be stated as ^ 6(j) 19|. It is well know that in the warm inflation scenario, the scalar 

held fluctuations are predominantly thermal rather than quantum ^ 201. In particular, 
for the weak dissipation regime, the amplitude of the scalar held fluctuation 50^ is given by 
6(y c:=L HT [20|. In this way, from Eqs.([8]), (ITT]) and (IT^ . the power spectrum Vh) results 


'P'R. — 


3\/^k / 


C. 


\QkC^ 


3-m 

l-m ll-3m , . , (3-m) / 2H^\ 2(4-m) 

4-m ( —fjj 4-m j \ -|- 

KT 


(19) 


Now, combining Eqs.([I]) and ffTT|) . we And that the power spectrum as function of the 
held (p can be written as 


1 _„ ^ 2/(4-m)+m-5 

= ki(p^ {B-^[K (p]) 


1 + 


2A‘^P 


kt{B ^[K (p])A'^ A 


3—m 

2(4-m) 


( 20 ) 


where the constant ki is defined as ki = 


Cs \ 


6kCo 


, ^ Q — zm , ^ m — o 

(Af) 4— (1 - /)4-m. 


By other hand, the scalar power spectrum can be expressed in terms of the number of 
e-folds N as 


1 — m 2/(4 —m) + m —5 

Vn{N) = k2 {B[J{N)]) ^ (J[iV]) 


1 + 


12712/^ 


kt{J[N]YA-A 


3—m 

2(4-m) 


( 21 ) 


where the quantities J{N) and k 2 are given by J{N) = 
respectively. 


l+fjN-l) 

Af 


, and k 2 = kiK 4 - 


From the definition of the scalar spectral index n^, given by the relation - 
and combining Eqs. (na, and fl20|) . the scalar spectral index can be written as 


1 _ d\n Vn 
^ dink ’ 


n.o = 1 


5 — m — 2/(4 — m) 
Af{A-m){B-^[K(P]y 


+ n2 + n2„ 


( 22 ) 
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where n 2 and are defined as 


and 


na = 


^3 = 


1-m 2{l - f){B-\K(t)\)-^/^ 


m 


SnAf 


1 + 




2A/(l-/)(3-m)^ 




-( 2 -/) 


Kr(S-i[i^0])2(W) 

2A‘^P 


-1/4 


1 + 


kt{B /) 


-1 


Kr(4 — m) 

Analogonsly as the scalar power spectrnm, the scalar spectral index can be expressed in 
terms of the nnmber of e-folds. Considering Eqs. flT7|l and ffTSjl . yields 

5 — m — 2/(4 — m) 


T) =1 — 

^ (4-m)[l + /(iV-l)] 


where now 


n2 = K 


1-m / 2(1 - /) {J[N])-f/^ 

4-mV nAf B[J{N)] 


1 + 


and 




1 + 


+ 712+ ns, 

\2A^P 

\2A^P 


(23) 


- 1/4 


-\ -1 


/ir(J[Ar])2(i-/) 


Kr(4 — m) 

It is well know that the generation of tensor pertnrbations dnring the inflationary scenario 
would produce gravitational waves. The power spectrum of the tensor perturbations is more 


complicated in our moc 


el because in brane-world gravitons propagate into the bulk. In this 


form, following Ref. 621], the power spectrum of the tensor perturbations Vg is given by 


Vg = 24:k{H/2'kYF‘^P). Here, x is defined as a: = Hrrip ■\/3/(47rr), and the function F{x) 
is given by 

F{x) = [\/l + a;2 — sinh“^(l/x)]“^'^^. 

The function F{x) is the correction to the standard General Relativity and arises from the 


normalization of a graviton zero-mode 


m\. 


In this way, we may compute an important observational quantity, the tensor-to-scalar 
ratio r, defined as r = Vg/Vn- Then considering Eq. fl20ll and Vg in terms of the scalar field, 
the tensor-to-scalar ratio can be written as 

QkA^P i-m , , . r 2A^p 


r(0) = 


vr^ ki 


0-4^ {B-\K4>\) 


3—m 

4 —m _j_ 


Kt{B ^[iC(/])2(l /) 


3—m 

2(4 —m) 


F\p. (24) 


Similarly, as before, the tensor-to-scalar ratio can be expressed in terms of the number 
of e-folds, yielding 

12A^P 


r(N) = P^(B[J(N)])-i^ (J|iV])- 


1 + 


Kr(J[iV]) 2 (i-/)_ 


3 — 771 
’ 2(4-m) 


F\N). (25) 
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FIG. 1: The evolution of the ratio T/2>H versus the scalar spectrum index Ug (upper panel) 
and the evolution of the ratio T/H versus the scalar spectral index Ug (lower panel) in the weak 
dissipative regime for the special case m = 1 (F oc T), for different values of the parameter 
In both panels, the solid, dotted, and dashed lines correspond to the pairs {A = 0.19, / = 0.31), 
{A = 0.28, / = 0.32) and(^ = 0.32, / = 0.30), respectively. Also, in these plots we have used the 
values = 70, nip = 1 and r = 10“^^. 

In Fig.(IT]) we show the evolution of the ratio F/3if versus the scalar spectral index (upper 
panel) and the evolution of the ratio T/H versus the the scalar spectral index (lower panel), 
in the weak dissipative regime for the specihc case in which the dissipative coefficient is given 
by F(0,T) = CfpT, i.e., m = 1. In both panels, we have used three different values of the 
parameter C^. The upper panel shows the condition for the weak dissipative regime in which 
F < 3H. In the lower panel we show the necessary condition for warm inflation scenario, 
in which the temperature T > H. Combining Eqs.([T]), 0131) and 014p . we numerically hnd 
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FIG. 2: The upper and lower panels show the evolution of the tensor-to-scalar ratio r versus the 
scalar spectral index Ug in the weak dissipative regime for the special case m = 0 (F oc </>), for 
different values of the parameter In both panels, the solid, dotted, and dashed lines correspond 
to the pairs (^4 = 0.47, / = 0.29), (^4 = 0.43, / = 0.28) and(^ = 0.39, / = 0.29), respectively. 
Also, in both panels we have used the values = 70, nip = 1, r = 10“^^, and the two-dimensional 
marginalized constraints from Ref. 


n 


50] (upper panel) and the new data from Planck 2015 M(i ower 


panel). 


the ratio F/Sif as a function of the scalar spectral index Ug. Analogously, from Eqs.ffT]) 
and dll. we numerically obtain the ratio between the temperature of the thermal bath T 
and the Hubble parameter i7, i.e., T/H m terms of the spectral index Ug. In both plots, 
we consider the values = 70, nip = 1 and r = 10“^. Also, numerically from Eqs. fl2T]) 
and (|2^ . we obtain the values A = 0.19 and / = 0.31 corresponding to the value of the 
parameter = 10“^, for the values V-ji = 2.43 x 10“®, Ug = 0.96, and the number of e-folds 
N = 60. Analogously, for the value = 10“^, the values obtained for the parameters A 
and / are given by A = 0.28 and / = 0.32, respectively. Finally,for the value = 10“®, 
we obtain the values A = 0.32 and / = 0.30. From the upper plot, we hnd an upper bound 
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for the parameter Ctji, from the condition of the weak dissipative regime, i.e., F < 3H, 
for which < 10“^. From the lower panel we obtain a lower bound for C^, considering 
the essential condition for warm inflation T > H, where > 10“®. In relation to the 
consistency relation r = r{ns), we hnd that the ratio r ~ 0 for this range of then 
the case m = 1 (or equivalently F oc T) during the weak dissipative regime is disproved 
from BICEP2 experiment, because r = 0.21 qq 5 and further the ratio r = 0 is discarded 
at 7.Oct. However, the Planck data analysis obtained only an upper limit for the tensor- 
to-scalar ratio r, given by r < 0.11, and then the range of is well corroborated from 
Planck satellite. In this form, for the specihc case of m = 1, the range of the parameter 

is given by 10“® < < 10“^. Also, we note that when we decrease the value of the 

parameter < 10“®, the value of the tensor to scalar ratio r ~ 0. In particular, for the 
value Cif) = 10“®, we get r |ns=o. 96 — 7 x lO”"^. It is interesting to note that the range for the 
parameter for the case F oc T results from the conditions F < 3H and T > H. 

In Fig.(j2]) we show the tensor-to-scalar ratio r versus the scalar spectral index Ug, for 
the special case of m = 0, i.e., F oc 0 in the weak dissipative regime (F < 3H). In the 
upper panel we show the two-dimensional marginalized constraints for the tensor-to-scalar 
ratio r = r{ns) (at 68% and 95% levels of conhdence), from the BICEP2 experiment data 
in connection with Planck satellite-f WP-|- highL50|. In the lower panel the new results 
from Planck 2015 [^|. Here, the marginalized joint 68% and 95% CL regions for the spectral 
index and ro,oo 2 - From Eqs. fl2^ and fl25l) . we numerically obtain the consistency relation 
r = r(n<j) and, as before, we consider three values of the parameter Cfp. Again, we take the 
values C-y = 70, nip = 1 and r = 10“^“^. From these plots we observe that the tensor-to 
-scalar ratio r ~ 0 for the specihc case of m = 0, and then this case is disproved from 
BICEP2 experiment (upper panel), however is well corroborated from Planck data and in 
special with the new data (lower panel). We observe that the new results from Planck 
2015 place more substantial limits on the consistency relation r = r{ns) compared with the 
BICEP2 experiment. 

In particular for the value of = 10“^^ (solid line in the hgure), we get that r |ns=o. 96 — 
2 X 10“^. From the essential condition for warm inhation T > H, we hnd a lower limit for 
the parameter given by Cip > 10“^^ and, considering the condition of weak dissipative 
regime, where F < 3H, we obtain un upper bound for C^, and it corresponds to < 10“^® 
(hgure not shown). In this form, for the special case of m = 0 (or equivalently F oc 0), the 
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range for the parameter is given by 10“^^ < C'c/i < 10“^®. 

For the cases m = 3 and m = —1, we hnd that the tensor-to-scalar ratio r ~ 0, then 
these cases are disproved from BICEP2 data, nevertheless are well corroborated from Planck 
satellite and Planck 2015 data. Considering the essential condition for warm inflation T > H, 
we obtain a lower bound for the parameter C^j)] for the case m = 3 the bound is 10^^ < 
and for the case m = —1 it corresponds to 10“^“^ < C^. In particular, for the value = 
10^^ corresponds to (^) |ns=o. 96 — 1-86 when m = 3, and for Cfj, = 10“^^ corresponds to 
|ns=o. 96 — 1-57 for the case m = —1. From the condition of weak dissipative regime 
F < 3H, we hnd an upper bound for for the case m = 3 the upper bound is < 10^^, 
and for the specihc case m = —1, results < 10“^®. In particular, for the value = 
10^^ corresponds to (^) |„^=o. 96 — 0.44 when m = 3, and for C^j, = 10“^® corresponds to 
(^) Us=o. 96 — 0.18 for the case m = —1. In this form, from the conditions T > H and 
F < 3H, the ranges for the parameter are given by: 10^^ < < 10^^ for the case 

m = 3, and 10“^^ < < 10“^® for the case m = —1. 


B. The strong dissipative regime. 


Now we consider the case of strong dissipative regime R> 1 (or equivalently F > 3H), 
together with the scalar factor a{t) of intermediate inhation, see Eq.([I]). Considering Eqs. 
(IH]) and flT^ . we hnd the solution for the scalar held 4’{t). In particular, we must to analyze 
our solution for two separate cases, namely m = 3 and m ^ 3. For the special case of m = 3, 
the solution for the scalar held results 




B[t] 
K ' 


(26) 


as before, the value of cf){t = 0) = 0o is an integration constant and iP is a constant given by 


K = 




r 


1/2 


2Ka 


3/8 


4-3(/-2) 

2 ( 1 -/) 

yJ 




( 1 -/)^ 

Af 


and B\t] is a new function that is dehned as 

B[t] = B 


2(Af) 


K-T ^2(1-/). 4/+ 3 15 


16(1-/)’ 16 


( 27 ) 


and corresponds to the incomplete beta function ^ . 
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On the other hand, the solution of the scalar held for the special case m 7 ^ 3 yields 

Bm[t] 


ip{t) -ipo = 


K„ 


( 28 ) 


here the scalar held 0 (f) is redehned as (pit) = and as before (po is an integration 

constant, that can be assumed (po = 0. The quantity Km is a new constant given by 


K — 


4Q 


1/2 




m/8 


4-m(/-2) / I_j.x4+m' 

( 2 ( 1 -/) ( (I-J) ' 

V 2 1 I 2 I/ 


for the special case in which m 7 ^ 3 is dehned as 

KT 


1 / 8 ( 1 -/) 


, and the new function Bm[t] 


Bm[t] = B 


2 ( 1 -/) 4:f + m 12+ m 

[ 2{Afy ’16(1-/)’ 16 


From Eqs.([I]), fl2B]) and fl25]) . we hnd that the Hubble parameter H = H{(j)) yields 

Af 


HU) = 

(H-i[iFln0])W 

and for the special case of m 7 ^ 3 we have 

Af 


, form = 3, 


«(/■) = 


(29) 


(30) 


(31) 


By considering Eq. flT^ , the scalar potential under the slow-roll approximation for both 
values of m is given by 


V{(j)) = T y-1 + 
for the special case m = 3, and 

H( 0 )=r (-1 + 


1 + 


2A^f 


Kr(H-i[iFln0])2{i-/) 


1 / 2 ^ 


(32) 


1 + 


2A^f 


1 l/ 2 ^ 


(33) 


for the value of m 7 ^ 3. 

The dissipative coefficient T in terms of the scalar held, can be obtained combining Eqs. 
(USD, ( 1 ^ - and ( 125 ]) to give 


r( 0 ) = S(j)-\B-^[K\n(j)])-^^ 


1 + 


2H^ 


KT 


-3/8 


Afjl-f) 


2kC~i 


for the case m = 3, in which the constant d is dehned as d = 
of m 7 ^ 3 we hnd 

2 ^ 21-"*/8 


3/4 


(34) 


For the value 




!(.(.) = 5„/‘-’"(By|A'„43l)-1 


1 + 


KT 


(35) 
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where 6m is a constant and is given by 6m = 


Afjl-f) 

2kC^ 


- m/4 


Analogous to the case of the weak dissipative regime, the dimensionless slow-roll param¬ 


eter e is given hj e = —^ = 


1 -/ 


A/(B-i [/s'In </.])/ 


, for m = 3 and for the value of m 7^ 3 we find 


£ = 


1-/ 


Again, if a > 0, then the scalar held (j) > exp[-bi?[(ij^)^/'f]], for the spe- 


Af(B^[KmV:>]y ""S'"*""’ - - - xxxx. OX.CXXCXX xxx^xxx yx / ^ [\ Af 

cial case m = 3, and for the case m 7^ 3 the new scalar held satished sp > 

Analogously as before, the inhationary scenario begins {e = 1) when the scalar held takes 
the value 0i = exp[i5[(irf)^/-^]], for m = 3, and ifi = for the special case 

of m 7^ 3. 

The number of e-folds N in this regime can be write using Eqs.(IT]), fl26|l . and fl28|) . to give 

/•i2 

Hdt = A[{B-^[K\n(P2]y - iB-yK\n(Pi]y], (36) 


N = 


Hi 


for the case of m = 3, and for the special case of m 7^ 3 we have 

N = Al(B-'lk„<p,]y - 


(37) 


Analogous to the case of the weak dissipative regime, now we will analyze the cosmological 
perturbations in which R = r/3H > 1 (strong dissipative regime). For the strong dissipative 
regime, the huctuations ^0^ is given by 6(j)‘^ ~ see Ref. 19| , where kp corresponds to 
the wave-number and it is determined as kp = \/VH = Hy/SR > H. In this way, the 
power spectrum of the scalar perturbation in this regime from Eqs. ®. oil and m , can 
be written as 


Pn ^ 


H 2 T 2 T 


K 


.Cfct>^)HP\-Hy-^ 


1 + 


2H^ 


(3m—6) 
16 


(38) 


4,r2 0 [ KT 

Analogously as before, we can hnd the power spectrum Vp in terms of the scalar held 0 
for both values of the parameter m. In this form, considering Eqs. ([I]), (126|1 . (128|1 and Eq. 
yields 

1 -3/16 


= k{B-\K\Yi 


3(5/-6) 


-3 


1 + 


2(A/)2(R-i[i^In0])-2d-/) 


KT 


(39) 


k 


for the case of m = 3. Here the constant k is given hj k = 


2kC-, 


11/8 


{AfY^/yi 


Now, the spectrum of the scalar perturbation for the special case of m 7^ 3, results 

2(4/)7B-‘l/?„.y>|)-"l'-7)' 


. , ~ 1 r ,, 3[/(2+m)-2m] 3(l-m) 

Vn = km{B;y[KY) -^-0^ 


1 + 


KT 


(3m—6) 
I6 


(40) 
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3771 + 2 


where km is a constant given by km = (2;^) * “ /)* • 

Also, the scalar power spectrum can be rewritten in terms of the number of e-folds N. 
From Eqs. fl36|l and ([37]), the power spectrum V-r becomes 


= k(jm ^exp ([1 + 


-3/16 


for the specihc case of m = 3, and for the case m 7^ 3 we have 


, r ,, 3[/(2+m)-2m] ~ , ^ ,,, 3(l-m) 

Vr = 7m{jm -^- {Bm[J[N]])^^ 


1 + 


2{AfyiJ[N) 


- 2 ( 1 -/) 


KT 


{3m—6) 
16 


(41) 


(42) 


3(1 —m) 

2Krr ' 


where the constant 7^ in the above equation is given by jm = km y 

For this regime the scalar spectral index Ug for the specihc case of m = 3 from Eqs. 
and (HOI) results 

ns = 1+ j {B~^[K\n(j)])~^ + rii + 712. 


8Af 


(43) 


1/2 

Here the new functions ni and n 2 are given by ni = —3 (|) ^ 


X (1 - /)V8 h 2(A/)^(h-Mh-in/.])-^(i-/) 


-2(1-/) 1 -1/16 


and 7/2 = 


0 

_ 3Af(l-f) I 5-1 




4«:t 


(5-1 [a: In 0])h 


-2 


X 




-1 


, respectively. By other hand, the scalar spectral index for 


the case m ^ 3 becomes 


3[/(7n + 2) — 2777] ~ f 

ns = l + - - ^ - -{Bm [Km^]) ^ + 772 „, 

where now the new functions Ui^ and 7/2^ are dehned by rii^ = (2)^ 


[4+m(/-2)] m-3 


X ^{Af)^{l - f)^{Bm^[Km^])-^^^<P^ 


2Af 

Z/D-lr 


^ 2{Aff{B-^[Km<p]) 


- 2 ( 1 -/) 


7 ^ 

771 — 4 


(44) 

—m/S 

and 


^2^ = 


(3m—6) 


4k;t 


(ll/)(l-/)(B 7 |/?„V>l)^ 


-2 


1 + 


2iAffiB-^lKm'fi]) 


- 2 ( 1 -/) 


-1 


Analogously as before, from Eqs. (l36|) and (j37|) . the scalar spectral index Ug can be rewrit¬ 
ten in terms of the number of e-folds as 


= ^ + ni + 772, 

0X1/ 


(45) 


for the specihc case of 777 = 3. Here the quantities 77i(J[iV]) and 772(J[iV]) are de- 

,1/2 1 / 1 X . .M/sf, . 2(lin2(,/(Vl'|-2(l-/)l"^/^® 


hned as 77i(J[iV]) = —3 (^) 


c-y2 y2KC-, 


(A/)-3/8(1-/)V8 


1 + 


2(Af)Ajm- 
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and n,(JlN]) = j[]V])/-2 + 

scalar spectral index for the value of m 7^ 3 yields 

3[/(m + 2) — 2m] 




-1 


, respectively. Analogously, the 


Us = 1 + 


(J[A^]) ^ + ni^+n2^, 


8Af 

[/(16-m.) + 2(m-6)] 


(46) 


[4+m(2-/)] 


where n,..(J|Af|) = (f) ^ (41/)-^^^^ 

[4+m(/-2)] , - ... ^ (3m-6) 


1 + 




X (J[]V])-^nwr^(B„[J|jV]|)-‘ and n2„(J|Af|) = 
2(a/)»(j[wl)-ai-ri 


(Am - f)(,JiN]y 


-2 


X 


-1 


Analogous to the case of the weak dissipative regime, the tensor-to-scalar-ratio r for the 
specihc case m = 3 can be expressed in terms of the scalar held to give 


r = -^(A/)2(rB-i[rAln0])T03 


1 + 


2(A/)2(5-i[.Aln0])-2d-/) 


KT 


16 


F2(0), (47) 


and the tensor-to-scalar-ratio for the case m ^ 3 yields 

6k 


r = 






1 + 




KT 


(3m—6) 
16 


F^iA- 


(48) 


Also, we obtain the tensor-to-scalar ratio in terms of the number of e-folds. In this form, 
combining Eqs. fl3^ and fH7|) we get 


r ^ ■!^(AfyiJlN])^ 

TT^K 


B[J[N]] 


K 


1 + 


2{Af)\J[N]) 


- 2 ( 1 -/) 


KT 


F^{N), (49) 


for the specihc case of m = 3. From Eqs. fl37p and the tensor to scalar ratio becomes 

^ 3(m—1) 

r = -^(A/)2(J[Ar])|[6"^+Aio-3™)-i6] ( ^-rn B^[J[N]] ^ ^ 




Km 


1 + 


2(rl/)7J|]V|) 


- 2 ( 1 -/) 


KT 


(3m—6) 
16 


F^N), 


(50) 


for the case of m 7^ 3. 

In Fig. on]) we show the evolution of the tensor-to-scalar ratio r and the ratio T/3H on the 
spectral index Ug in the strong dissipative regime, for the special case in which E = Cy,T^/(jA, 
i.e., m = 3. In both panels we have used three diherent values of the parameter C^). In 
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FIG. 3: The evolution of the tensor-to-scalar ratio r (upper panel) and the evolution of the 
ratio T /‘iH versus the scalar spectral index Ug (lower panel) in the strong dissipative regime, for 
the special case m = 3 (F oc T^/cfP') for three different values of the parameter C^. In both 
panels, the solid, dotted, and dashed lines correspond to the pairs {A = 1.32 x 10“®, / = 0.75), 
{A = 7.68 X 10“®, / = 0.87) and {A = 5.40 x 10“®, / = 0.97). In both panels we have used 
the values = 70, rup = 1, r = 10“^^, and the two-dimensional marginalized constraints from 
Ref. 50]. 


the upper panel, as before, we show the two-dimensional constraints in the n^-r plane from 
BICEP2 and Planck data (68% and 95% levels of conhdence) [^ . In the lower panel we 
show the evolution of the ratio P /3H during the warm inflation scenario, and we corroborate 
that our model satishes the strong dissipative regime, i.e., P > 3H. In order to write down 
the consistency relation r = r{ns), we numerically hnd from Eqs. (14511 and fl49ll . the relation 
r = rijis) (upper panel). Similarly, from Eqs. fIM]) . fimi and (145]) we obtain the ratio V/3H 
as a function of the spectral index rig (lower panel). In these plots we consider the values of 
= 70, nip = 1 and r = Also, numerically from Eqs. (l4TD and (145]) for the special 

case of m = 3, we obtain the values A = 1.32 x 10“^ and / = 0.75 corresponding to the value 
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= 10®, for P-ji = 2.43 x 10“®, Ug = 0.96 and N = 60. Analogously, for the parameter 
6*0 = 10®, we numerically find the the values A = 7.68 x 10“® and / = 0.87. Finally, for 
the value = 10^ we obtain A = 5.40 x 10“® and / = 0.97. From the upper panel we 
obtain that the range for the parameter is 10^ < < 10®, which is well supported from 

observational data. However, from the lower panel we observe that for values of < 10® 
the model is disapproved from the condition the strong dissipative regime, since the ratio 
F/3i7 < 1. Here, we note that from the condition of the strong dissipative regime, i.e., 
F > 377, we have obtained a lower bound for the parameter C^. In this form, for the value 
m = 3 ( in which F oc /(fP') the range for the parameter is given by 10® < < 10®, 

which is well supported from observational data together with the conditions of the strong 
dissipative regime F > 3H and T > H. 

During the strong dissipative regime for the specihc case of m = 1 (F oc T) , we obtain 
that for the value of the parameter > 10“^ the model is well corroborated from the 
condition F > 3H and the necessary condition for warm inflation T > H (hgure not shown). 
For the tensor-to-scalar ratio, we hnd that the ratio r ~ 0 for this lower bound, then the 
case m = 1 is disproved from BICEP2 experiment, since the ratio r = 0 is discarded at 7. Oct. 
However, from the Planck data, the value > 10“^ is well supported. Also, we note that 
when we increase the value of the parameter > 10“^, the value of the tenso-to-scalar 
ratio r ~ 0. 

For the cases m = 0 (F oc 0) and m = — 1 (F oc (jP/T), we hnd that these models in the 
strong dissipative models are disproved from BICEP2 and Planck data, because the scalar 
spectral index > 1, and then these models do not work. 


III. CONCLUSIONS 

In this paper we have studied warm-intermediate inhation in the context of Randall- 
Sundrum H brane-world cosmological model. Considering the slow-roll approximation dur¬ 
ing the weak and strong regime, we have obtained analytical solutions of the full effective 
Friedmann equation for a hat Universe in this brane-world model. Here we have consid¬ 
ered a standard scalar held 0 together with a general form of the dissipative coefficient 
F oc T™/0™“^. In special, we analyzed the values m = 3, m = l, m = 0 and m = —1, 
that can be have found in the literature for this dissipation coefficient. Studying the weak 
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and strong dissipative regimes, we have obtained analytical expressions for the appropriate 
Hnbble parameter, effective potential, scalar power spectrnm, scalar spectral index and the 
tensor-to-scalar ratio. Dnring both regimes we have studied the slow roll analysis and we 
compared with the two-dimensional marginalized constraints (68% and 95% C.L.) r = r{ns) 
plane from observational data. Also, we have obtained a constraint for the parameter 
(see Eqj6]) from BICEP2 and Planck 2015 data together with the essential condition for 
warm inflation T > H and the condition from the weak P < 3H (or strong P > 3H) regime. 

For all the models (different values of the parameter m) in the weak dissipative regime, 
we have found a lower bound for the parameter from the essential condition for warm 
inflation, in which the temperature of the thermal bath T > H. Also, we have obtained 
an upper bound for from the condition P < 3H, i.e., the weak dissipative regime. 
Additionally, we have observed that the consistency relation r = r{ns) ~ 0, in the weak 
dissipative scenario, and the models are disproved from BICEP2, but are well corroborated 
from Planck satellite, since r < 0.11. 

For the strong dissipative scenario, we have found that the range for the parameter is 
given by 10® < < 10® in the specihc case of m = 3 i.e., P oc Here, we have found 

an upper bound from BICEP2-Planck r — Ug plane and a lower bound from the condition 
of the dissipative regime in which P > 3H, also in this range of the necessary condition 
for warm inflation T > if is satished. For the case m = 1 (P oc T), we have obtained that 
the tensor-to-scalar ratio r ~ 0, and also we have found a lower limit for the parameter 
from the condition P > 3H. Finally, we have found that for the cases m = 0 and m = —1, 
these warm-intermediate inflationary models are disproved from observational data, since 
the scalar spectral index Ug > 1, then the models P oc 0 and P oc 0^/T do not work. 

Acknowledgment s 

R.H. was supported by Comision Nacional de Ciencias y Tecnologia of Chile through 
FONDECYT Grant N® 1130628 and DI-PUCV N® 123.724. N.V. was supported by Comision 
Nacional de Ciencias y Tecnologia of Chile through FONDECYT Grant N® 3150490. 


[1] A. Guth , Phys. Rev. D 23, 347 (1981) 


22 



[2] A.A. Starobinsky, Phys. Lett. B 91, 99 (1980) 

[3] A.D. Linde, Phys. Lett. B 108, 389 (1982) 

[4] A.D. Linde, Phys. Lett. B 129, 177 (1983) 

[5] A. Albrecht and P. J. Steinhardt, Phys. Rev. Lett. 48,1220 (1982) 

[6] K. Sato, Mon. Not. Roy. Astron. Soc. 195, 467 (1981). 

[7] V.F. Mukhanov and G.V. Chibisov , JETP Letters 33, 532(1981) 

[8] S. W. Hawking,Phys. Lett. B 115, 295 (1982) 

[9] A. Guth and S.-Y. Pi, Phys. Rev. Lett. 49, 1110 (1982) 

[10] A. A. Starobinsky, Phys. Lett. B 117, 175 (1982) 

[11] J.M. Bardeen, P.J. Steinhardt and M.S. Turner, Phys. Rev.D 28, 679 (1983). 

[12] D. Larson et ai, Astrophys. J. Suppl. 192, 16 (2011). 

[13] C. L. Bennett et ai, Astrophys. J. Suppl. 192, 17 (2011) 

[14] N. Jarosik et ai, Astrophys. J. Suppl. 192, 14 (2011) 

[15] G. Hinshaw et al. [WMAP Gollaboration], Astrophys. J. Suppl. 208, 19 (2013) 

[16] P. A. R. Ade et al. [Planck Gollaboration], Astron. Astrophys. 571, A16 (2014) 

[17] P. A. R. Ade et al. [Planck Gollaboration], Astron. Astrophys. 571, A22 (2014) 

[18] P. A. R. Ade et al. [Planck Gollaboration], arXiv:1502.02114 [astro-ph.CO]. 

[19] A. Berera, Phys. Rev. Lett. 75, 3218 (1995); A. Berera, Phys. Rev. D 55, 3346 (1997). 

[20] L.M.H. Hall, LG. Moss and A. Berera, Phys.Rev.D 69, 083525 (2004). 

[21] A. Berera, Phys. Rev.D 54, 2519 (1996). 

[22] A. Berera, 1. G. Moss and R. O. Ramos, Rept. Prog. Phys. 72, 026901 (2009); M. Bastero-Gil 
and A. Berera, Int. J. Mod. Phys. A 24, 2207 (2009). 

[23] A. Sen, JHEP 0204, 048 (2002). 

[24] K. Akama, Lect. Notes Phys. 176, 267 (1982); V. A. Rubakov and M. E. Shaposhnikov, Phys. 
Lett. B 159, 22 (1985); N. Arkani Hamed, S. Dimopoulos, and G. Dvali, Phys. Lett. B 429, 
263 (1998); M. Gogberashvili, Europhys. Lett. 49, 396 (2000); L. Randall and R. Sundrum, 
Phys. Rev. Lett. 83, 3370 (1999); 83, 4690 (1999). 

[25] T. Shiromizu, K. Maeda, and M. Sasaki, Phys. Rev. D 62, 024012 (2000). 

[26] P. Binetruy, C. Deffayet, and D. Langlois, Nucl. Phys. B565, 269 (2000); P. Binetruy, C. 
Deffayet, U. Ellwanger, and D. Langlois, Phys. Lett. B 477, 285 (2000). 

[27] L. Randall and R. Sundrum, Phys.Rev.Lett.83, 4690 (1999). 


23 


[28] L. Randall and R. Sundrum. Phys. Rev. Lett., 83, 3370, (1999). 

[29] W. D. Goldberger and M. B. Wise, Phys. Rev. Lett. 83, 4922 (1999); W. D. Goldberger and 
M. B. Wise, Phys. Lett. B 475, 275 (2000. 

[30] O. DeWolfe, D. Z. Freedman, S. S. Gubser and A. Karch, Phys. Rev. D 62, 046008 (2000); 
W. D. Goldberger and 1. Z. Rothstein, Phys. Lett. B 491, 339 (2000); N. Arkani-Hamed, 
S. Dimopoulos and J. March-Russell, Phys. Rev. D 63, 064020 (2001). M. A. Luty and R. 
Sundrum, Phys. Rev. D 64, 065012 (2001); R. Hofmann, P. Kanti and M. Pospelov, Phys. 
Rev. D 63, 124020 (2001); J. Garriga, O. Pujolas and T. Tanaka, Nucl. Phys. B 605, 192 
(2001); A. Flachi, 1. G. Moss and D. J. Toms, Phys. Rev. D 64, 105029 (2001); S. Nojiri, S. 
D. Odintsov and S. Zerbini, Glass. Quant. Grav. 17, 4855 (2000); 1. Brevik, K. A. Milton, S. 
Nojiri and S. D. Odintsov, Nucl. Phys. B 599, 305 (2001). 

[31] T. Tanaka and X. Montes, Nucl. Phys. B 582, 259 (2000); S. Mukohyama and L. Kofman, 
Phys. Rev. D 65, 124025 (2002). 

[32] R. Maartens, D. Wands, B. A. Bassett, and 1. P. C. Heard, Phys. Rev. D 62, 041301 (2000). 

[33] J. M. Gline, C. Grojean, and G. Servant, Phys. Rev. Lett.83, 4245 (1999). 

[34] P. Brax and G. van de Bruck, Glass. Quant. Grav. 20, R201 (2003); T. Clifton, P. G. Ferreira, 
A. Padilla and C. Skordis, Phys. Rept. 513, 1 (2012). 

[35] C. Csaki, M. Graesser, C. F. Kolda and J. Terning, Phys. Lett. B 462 34 (1999); D. Ida, 
JHEP 0009, 014 (2000); R. N. Mohapatra, A. Perez-Lorenzana, and C. A. de S. Pires, Phys. 
Rev. D 62, 105030 (2000); R. N. Mohapatra, A. Perez-Lorenzana, and C. A. de S. Pires, Int. 
J. Mod. Phys. A 16, 1431 (2001); Y. Gong, arXiv:gr-qc/0005075; R. Herrera, Phys. Lett. B 
664, 149 (2008); Y. Ling and J. P. Wu, JCAP 1008, 017 (2010). 

[36] R. Maartens, arXiv:gr-qc/0101059. 

[37] F. Lucchin and S. Matarrese, Phys. Rev. D32, 1316 (1985). 

[38] J. D Barrow, Phys. Lett. B 235, 40 (1990); J. D Barrow and P. Saich, Phys. Lett. B 249, 406 
(1990);A. Muslimov, Class. Quantum Grav. 7, 231 (1990); A. D. Rendall, Class. Quantum 
Grav. 22, 1655 (2005). 

[39] J. D Barrow and A. R. Liddle, Phys. Rev. D 47, R5219 (1993); A. A. Starobinsky JETP 
Lett. 82, 169 (2005); S. del Campo, R. Herrera, J. Saavedra, C. Campuzano and E. Rojas, 
Phys. Rev. D 80, 123531 (2009); R. Herrera and E. San Martin, Eur. Phys. J. C 71, 1701 
(2011); R. Herrera and M. Olivares, Mod. Phys. Lett. A 27, 1250101 (2012); R. Herrera and 


24 



M. Olivares, Int. J. Mod. Phys. D 21, 1250047 (2012). 

[40] W. H. Kinney, E. W. Kolb, A. Melchiorri and A. Riotto, Phys. Rev. D 74, 023502 (2006); R. 
Herrera and E. San Martin, Int. J. Mod. Phys. D 22, 1350008 (2013). 

[41] J. D. Barrow, A. R. Liddle and C. Pahud, Phys. Rev. D, 74, 127305 (2006); R. Herrera, M. 
Olivares and N. Videla, Eur. Phys. J. C 73, 2295 (2013). 

[42] T. Kolvisto and D. Mota, Phys. Lett. B 644, 104 (2007); Phys. Rev. D. 75, 023518 (2007). 

[43] I. Antoniadis, J. Rizos and K. Tamvakis, Nucl.Phys. B 415, 497 (1994). 

[44] D. G. Boulware and S. Deser, Phys.Rev. Lett. 55, 2656 (1985); Phys. Lett. B 175, 409 (1986). 

[45] S. Mignemi and N. R. Steward, Phys. Rev. D 47, 5259 (1993); P. Kanti, N. E. Mavromatos, 
J. Rizos, K. Tamvakis and E. Winstanley, Phys. Rev. D 54, 5049 (1996); Ch.M Chen, D. V. 
Gal’tsov and D. G. Orlov, Phys. Rev. D 75, 084030 (2007). 

[46] S. Nojiri, S. D. Odintsov and M. Sasaki, Phys. Rev. D 71, 123509 (2004); G. Gognola, E. 
Eizalde, S. Nojiri, S. D. Odintsov and E. Winstanley, Phys. Rev. D 73, 084007 (2006). 

[47] A. K. Sanyal, Phys. Lett. B, 645,1 (2007). 

[48] 1. Antoniadis, J. Rizos and K. Tamvakis, Nucl. Phys. B. 415, 497 (1994); S. Nojiri, S. D. 
Odintsov and M. Sasaki, Phys. Rev. D. 71, 123509 (2004); G. Gognola, E. Eizalde, S. Nojiri, 
S. D. Odintsov and E. Winstanley, Phys. Rev. D. 73, 084007 (2006). 

[49] R. Herrera, N. Videla and M. Olivares, Phys. Rev. D 90, no. 10, 103502 (2014); R. Herrera, 
M. Olivares and N. Videla, Int. J. Mod. Phys. D 23, no. 10, 1450080 (2014); R. Herrera, 
M. Olivares and N. Videla, Phys. Rev. D 88, 063535 (2013); R. Herrera and E. San Martin, 
Int. J. Mod. Phys. D 22, 1350008 (2013). R. Herrera and E. San Martin, Eur. Phys. J. C 71, 
1701 (2011). 

[50] P. A. R. Ade et al. (BICEP2 Collaboration), Phys. Rev. Lett. 112, 241101 (2014); Astrophys. 
J. 792, 62 (2014). 

[51] R. Adam et al. [Planck Collaboration], arXiv: 1409.5738 [astro-ph.CO]. 

[52] M. A. Cid, S. del Campo and R. Herrera, JCAP 0710, 005 (2007). 

[53] 1. G. Moss and C. Xiong, arXiv:hep-ph/0603266. 

[54] A. Berera, M. Gleiser and R. O. Ramos, Phys. Rev. D 58 123508 (1998). 

[55] A. Berera and R. O. Ramos, Phys. Rev. D 63, 103509 (2001). 

[56] Y. Zhang, JCAP 0903, 023 (2009). 

[57] M. Bastero-Gil, A. Berera and R. O. Ramos, JCAP 1107, 030 (2011). 


25 



[58] M. Bastero-Gil, A. Berera, R. O. Ramos and J. G. Rosa, JCAP 1301 , 016 (2013); M. Bastero- 
Gil, A. Berera, R. O. Ramos and J. G. Rosa, arXiv: 1404.4976 [astro-ph.GOj. 

[59] S. Bartrum, M. Bastero-Gil, A. Berera, R. Cerezo, R. O. Ramos and J. G. Rosa, Phys. Lett. 
B 732 , 116 (2014). 

[60] J. Yokoyama and A. Linde, Phys. Rev D 60, 083509, (1999); R. Herrera, M. Olivares and 
N. Videla, Phys. Rev. D 88, 063535 (2013). 

[61] M. Abramowitz and 1. A. Stegun, Handbook of Mathemati- 733 cal Functions with Formulas, 
Graphs, and Mathematical 734 Tables, 9th ed. (Dover, New York, 1972); P. Appell and 735 
J. Kampe de Feriet, Fonctions Hypergeeomeetriques et 736 Hyperspheeriques: Polynomes 
DHermite (Gauthier-Villars, 737 Paris, 1926); A. Erdelyi, Acta Math. 83 , 131 (1964); H. 
738 Exton, in Multiple Hypergeometric Functions and Appli- 739 cations (Wiley, New York, 
1976). 

[62] D. Langlois, R. Maartens and D. Wands, Phys. Lett. B 489 , 259 (2000). 


26 


